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Abstract 



O-i, We use a path integral formalisni to derive the semiclassical series for the 

■y \ partition function of a particle in D dimensions. We analyze in particular 

C^ ' the case of attractive central potentials, obtaining explicit expressions for the 



fluctuation determinant and for the semiclassical two-point function in the 
special cases of the harmonic and singe-well quartic anharmonic oscillators. 
The specific heat of the latter is compared to precise WKB estimates. We 
rS \ conclude by discussing the possible extension of our results to field theories. 
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I. INTRODUCTION 

As is well known [0-|^ , the partition function of a particle of mass m interacting with a 
potential V{x.) and a thermal reservoir at temperature T can be written as a path integral 

{(3 = i/kBTy. 

Z{(3)= f t/^xop(/?;xo,xo), (la) 

JlRO 



p(/?;xo,Xo)= / [I?x(r)] e-^M/\ (lb) 



S\x] = / dr 
Jo 




(Ic) 



This path integral may be approximated in a number of ways: depending on the circum- 
stances, one may resort to perturbation theory around the exactly soluble harmonic oscil- 
lator, variational estimates, or lattice Monte Carlo calculations (such techniques carry over 
to Quantum Statistical Field Theory, where free fields play the role of unperturbed uncou- 
pled harmonic oscillators). Semiclassical techniques can also be used in approximating this 
integral. It is their virtues and shortcomings in applications to Statistical Mechanics that 
we intend to discuss. 

Semiclassical techniques have proven extremely important in the discussion of the tran- 
sition from Quantum to Classical Mechanics PJ^. In the present context, however, we shall 
use them in the opposite sense: to systematically incorporate fluctuations (thermal and 
quantum) to a description that has one or more solutions of the "Euclidean" equations of 
motion as its starting point. (Heretofore, we call these solutions "trajectories" or "classical 
paths".) The Euclidean character is of crucial importance: first, it restricts the trajectories 
to be global minima of the Euclidean action [^ — any others are exponentially suppressed; 
in addition, it leads to classical mechanics problems whose potencial is minus the physical 
one. Since we are interested in traces of operators, only closed trajectories will contribute. 
All this dramatically reduces the number of trajectories. In the specific examples of the 
harmonic and single-well quartic anharmonic oscillators, only one trajectory exists once the 
initial position and "time-of-fiight" f3h are fixed. 

Thanks to the features described in the previous paragraph, in a recent paper 0] we were 
able to construct the full semiclassical series for the partition function of a particle in one 
dimension from the mere knowledge of the trajectories. We obtained fluctuation determi- 
nants in a straightforward manner, by-passing the solution of the equivalent boundary- value 
problems, generated all the terms of the series in a systematic way, and could show that 
each term has a non-perturbative character, corresponding to sums over infinite subsets of 
perturbative graphs. Furthermore, we showed that the construction actually contains all 
the perturbative diagrams and many more. As an application of the method, we evaluated 
the ground-state energy and the specific heat of the single-well quartic anharmonic oscil- 



lator, achieving for the former good agreement with precise numerical results [0, and for 



the latter a result which has the correct high temperature limit, in contrast with the one 
obtained via conventional perturbation theory around the minima of the potential. 

In this article, we present the D-dimensional generalization of the method. Indeed, we 
are able to prove, as in the one-dimensional case, that it is possible to evaluate each term of 
the semiclassical series for the partition function using the classical path(s) as the only input. 
For the sake of simplicity, we concentrate on the case of attractive central potentials. In 
such potentials, as will be shown below, the only trajectories that contribute to the partition 
function are the ones with zero angular momentum. The discussion for arbitrary potentials 
is left for an appendix. As examples, we consider the isotropic harmonic oscillator and the 
single-well quartic anharmonic oscillator; in particular, we compute the specific heat of the 
latter in the lowest order semiclassical approximation for a few values of the temperature 
and for D = 1, 2 and 3. 

The article is organized as follows: Section II presents the derivation of the semiclassical 
series for a generic potential in an arbitrary number of dimensions, and the explicit formulae 
for the fluctuation determinant and the semiclassical two-point function in the particular 
case of attractive central potentials; Section III illustrates these results in the cases of the 
harmonic and single-well quartic anharmonic oscillators; Section IV presents our conclusions. 
In the Appendix, we show how to obtain the fluctuation determinant and the semiclassical 
two-point function in the case of an arbitrary potential in D dimensions. 

II. THE SEMICLASSICAL EXPANSION IN STATISTICAL MECHANICS 

A. General formalism 

The procedure to generate a semiclassical series for Z{(3), Eq. (|1|), was carried out in 
detail in Ref. |^ for the one- dimensional case {D = 1). Here we shall only sketch its 
generalization for arbitrary D (for a detailed discussion of the semiclassical expansion in 
quantum mechanics using path integrals, see Refs. [0,0). The first step is to find the 
minima Xc(r) of the Euclidean action 5'[x]. They satisfy the Euler-Lagrange equation 



mxr 



V\/(xJ = , (2) 



subject to the boundary conditions Xc(0) = Xc(/3/i) = xq; for simplicity, we shall assume 
here that there is only one minimum. The next step is to functionally expand the Euclidean 
action around it. Writing x(r) = ^c{t) + u(r), with u(0) = u{ph) = 0, we have ^[x] = 
S[^c] + 5'2[u] -|- (JS'fu], where 



1 f/^f^ 

5*2 [u] = - / druiir) 

2 Jo 



2 



7 



Uj{t) , (3a) 



SS[u]= / dT6V{r,u)= / rfr^-a,,...a,„l^(x,)M,,(r)...n,„(r); (3b) 

Jo Jo ^^ n\ 

the indices i, j, . . . run from 1 to D, and repeated indices are summed. Inserting this de- 
composition of S into (^ and expanding e~^^l^ in a power series yields the semiclassical 
expansion of 2'(/5): 



z(,,./^„.^.e--/7::rw,,.-^.-.£i(-^)". (, 



The first term of tlie series corresponds to tlie quadratic approximation to tlie partition 
function, which we denote by Z2{P): 



JiRO Ju{0)=0 

= I d^xoe-^[^d/"A-l/^ (5) 

J]RD 

where A is the determinant of the fluctuation operator JF: 

(P 
A = Det J", J^ij = -m — 6ij + didjV{yic) ■ (6) 

The other terms of the series (|) lead to integrals of the type 

/•u(/3a)=0 , , , 

K(r0...n.,(r,))= / [I?u(r)] e-^^M/^n,,(ri) • • • n,,(r,) . (7) 

Ju(0)=0 

Since the action 52 [u] is quadratic, one can show that 

p 

if k is even, and zero otherwise. X^p denotes the sum over all possible pairings of the Tj^., 
and Gijij, r') is the solution of 



rf2 
-m -— 6ij + didjV{yic 



G',fe(r,r') = 5,fc5(r-r'), (9) 



satisfying the boundary conditions 

G',fc(0,r') = G',fe(/3;i,r') = 0. (10) 

In the Appendix we present a recipe for obtaining A and G^ (r, r') using the general solution 
of the equation of motion (^ as the only input. 

B. Central potentials 

To illustrate the formalism of the previous section, let us apply it to the case of central 
potentials, i.e., V = V{r), where r = |x|. First of all, we note that, because of the radial 
symmetry, p(/5; xq, xq) can only depend on tq = |xo|. Thus, without loss of generality, we 
may take xq = r^ei, where ei is the unit vector pointing in the Xi-direction, and perform 
the angular integration in ([Ta|) to obtain 

27r-^/2 roc 

^^^') = Y{DjY) L drort'pi(^;roe,;roe^). (11) 



In general, there are many classical trajectories satisfying the boundary conditions x(0) = 
x(/5^) = vq Si. However, they are all radial if the potential is purely attractive [i.e., V'{r) > 
for r > 0]. Indeed, in this case the Euclidean motion is equivalent to that of a particle in a 
repulsive central potential, so that a closed classical trajectory necessarily has zero angular 
momentum. Besides, this trajectory is unique if the potential is smooth at the origin, i.e., 
V'{0) = 0. 

For a trajectory lying in the Xi-axis, ^c{t) = rci'^)^!, the fluctuation operator JF is 
diagonal in the indices i,j. Indeed, since V = V{r), we have 



V'(r) 



V'ir)-^ 



Jb'iJb -i 



(12) 



which, for Xi = VcSn, gives didiV{rc) = V"{rc), didiV{rc) = r^^ V'{rc) ioi i = 2, . . . ,D, and 
didjV{rc) = OiU^j. Thus, A = A^ Af ~\ where 

Ae = Det [-md'^ + V"{r,)] , A^ = Det [-mdl + r^^ Vir^)] (13) 

(£ and t stand for longitudinal and transverse, respectively). 

The Green's function Gij also becomes diagonal in this case: Gu = Gi, Gu = Gt for 
i = 2, . . . ,D, and Gij = ii i ^ j, where 

[-m 9^ + V"{r,)] Gfir, r') = Sir - r') , (14a) 

[-m dl + r;^ V'{r,)] ^^(r, r') = 5{t - r') . (14b) 

A^ and Gi{t, t') are the fluctuation determinant and semiclassical Green's function that 
appear in the one-dimensional version of the problem, which was studied in detail in Ref. 
PJ. There, the following results were derived: 

A 2vr/i , l],(0,r<)l],(r>,/3;^) 

m mUi[0, pn) 

where T^{Ty) = min(max){r, r'} and 

i2^(r,r) = — — — -— —— — — - , (16) 

Va{T')Vb{r')~Va{r')Vb{r') 

with rjai^r) and rjb^r) any two linearly independent solutions of the homogeneus equation 

[-md^ + V"{r,)]r]{T)=0. (17) 

By differentiating the equation of motion mr^. — V'{r() = with respect to r, one can 
verify that rja^r) = fdr) is one such solution. The other can be taken as [^ Vbi.^) = 
Tcij) Jq dr' [rc(r')]^^. For such a choice, the denominator of ^^(r, r') is equal to 1 and, since 
r]b{0) = 0, one has A^ = {27ih/m) r]a{0)rih{Ph). Because of these simplifying features, we 
shall refer to those solutions as the "canonical" solutions of (0). 

In order to obtain At and Gt{r, r') one simply replaces VLi{t, t') in (^) by 



O / /N ^ 'Pair) M^') - fai^') M^) f.nS 

'''^^'^^- Mr') Mr') -Mr') Mr')' ^ ^ 

where (pa(r) and (pb{r) are two linearly independent solutions of 

l-m dl + r-^ V'{r,)\ M) = • (19) 

It immediately follows from the equation of motion that ipair) = r^r) is one such solution. 
Another one is ifb{r) = rc{T) Jq dr' [rc{T')]~'^. They form a pair of canonical solutions of 
(0)- 

III. APPLICATIONS 

Using the results of the previous section we may write the quadratic approximation to 
the partition function as 

T{D/2) 

It can be readily calculated from the knowledge of Xc(r) alone. This will be acomplished 
below for both the harmonic and single-well quartic anharmonic oscillators. 



^'^^^ = ffm r^r^rt'e-'^-^y^^ (A, Af-)"^/^ . (20) 



n D/2 »oo 

ZiP) = ^^ I dr, rt' e-^[-l/n a-^/^ . (22) 



A. The harmonic oscillator 

As a first example, we consider the D-dimensional (isotropic) harmonic oscillator, 

l/(r) = -mcjV. (21) 

Since the potential is quadratic, 6V{t, u) = and Z{P) = Z2{P). Besides, r~^ V'{r) = V"{r), 
so that At = Ae. Thus, 

27r^/^ 
T{D/2) Jo 

The solution of the equation of motion is straightforward, and yields 

_ rp cosh[^(r - ph/2)] 
'^^^^^" cosh{f3nu/2) • ^^^^ 

The classical action can be readily computed, giving 

S[rc] = mujrl taxi\i{(3huj / 2) . (24) 

As solutions of (P!7| ) we may take Tjair) = cosh(ci;r) and i]h{T) = sinh(ct;r). This gives 
Qi{T,T') = uj~^ sinh[ti;(r' — r)], so that 

A, = ^^^ ^i^M/^^c.) _ 

Inserting ( p4D and ( pSj) into (^) and performing the integral, we obtain 

Z{p) = [2 smh{(]huj/2)y^ , (26) 

which is the well-known result for the partition function of the D- dimensional harmonic 
oscillator. 



B. The single-well quartic anharmonic oscillator 



Let us now consider the potential 



V{r) = - mu^r^ + ^ Ar^ (A > 0) . 



(27) 



In order to simplify notation, it is convenient to replace r and r by g = (X/muj^Y''^ r and 
6 = LOT, respectively. In the new variables, the equation of motion reads 



(Pq o 



(28) 



whose solution, taking into account the boundary conditions, is 

q,{9)=qtnc{ue,k), (29) 

where nc(u, k) = l/cn(M, k) is one of the Jacobian Elliptic functions ||TB|-|T3|, and 



Ug 



l + q! 



e 



k 



2 + q! 



\2(l + g| 



2\ ' 



where = phu. The relation between go and qt is obtained by taking 6* = in (p9|): 



go = gc(0) = qt^cue . 

(From now on we shall omit the /c-dependence in the Jacobian Elliptic functions.) 
The classical action can be written as S[rc\ = {m'^uj^/X) I[qc], where 



(30) 



(31) 



i[q] = / de 

Jo 



r + U{q) 



U{q) = \q' + \q' 



Using \q1 — U{qc) = —U{qt), we may rewrite /[gj as 



nqc] = U{qt) + 2 /'" dq p [U{q) - U{qt)] 
J lit 



(32) 



(33) 



Performing the integration and using ( PT| ) yields 

I[qc] = Q{^^q't+\qi)+\{-^^Wt 



E{ipe,k) + -q^ue 



+ snwe (1 + -g^^nc^Mej Jl + -qf {1 + nc'^ue) 



(34) 



where E{ip, k) denotes the Elliptic Integral of the Second Kind and ipe = arccos[gc(6')/go] 
arccos(cnu0). 

The canonical solutions of (|17[) and (pJl) are given by 



■ snuednug 



r]a{e) = uqt \jl + g| 

1 sn-uediiMg 



ca^ue 



(35a) 



vm 



uj^qt (1 + qf) Clique 



en ug dn ug 



sn Ug 

MO)- "^"^ 



+ ik' - 1) 



sn ue en M51 



(^-0) 



u;/ 



k\ Jl + g| 



dnue 

E{^g,k) + {e-l)ug-{9^0) 



(35b) 

(35e) 
(35d) 



Thus, 



1 - P 2A;2 - 1 ^ 

Me H 7-^ — E(v90,/i;) 



fc2 



47r^ sn^Me dn Me 

"^^ y^l + g2 cn%e 

en Me dn Me , ,^, snweenwe 

H ^^ + (1 - F) r^ 

sn Me dn uq 

iiih ne^Me 



fc2 



"^^ k^Jl + qi 



E{ifQ,k) + {e-l)ue 



(36a) 
(36b) 



We now have all the neeessary ingredients to eompute the quadratic approximation to 

z{py. 



Z2{P) 



T{D/2) \ A 



'^y fdqoqre-^^^^y^{A,Ar' 



-1/2 



(37) 



where g = h\/m^uj^. However, to perform the integral over go one must write /[^c], A^ and 
Ai in terms of go- In view of Eq. (|31|), it is much simpler to change the variable of integration 
from go to g^: 



Z2{P) 



27r^/^ 
r(D/2) V A 



where ge = limg^^oo qt{,qa, ©)• The Jacobian {dqo/dqt)e can be obtained directly from (|3TD 
by differentiation or, more simply, by using the identity [§ 



muj 



U'{qt)At 



'dqo\ 

.^^Je ^vr/i /2 [f/(go) _ U{qt)] 



(39) 



As an application, we may use ( |38D to calculate the specific heat of the D-dimensional 
single-well quartic anharmonic oscillator, given by 



C = (3 



(92 



(40) 



Using the program MAPLE, we computed this expression for a few values of the temperature 
17|. In Fig. [l|, we present the results for D = 1, 2 and 3. In Fig. |^, we compare the 



semiclassical approximation with: i) the classical result, in which the partition function is 
given by 

ii) with the lowest order WKB approximation, in which the energy levels entering the ex- 
pression Z = J2n e~^^" are given by the Bohr- Sommerf eld formula 

^ ^2m [En - V{x)] dx=(^n+^^h (n = 0, 1, 2, . . .); (42) 

and iii) with the specific heat of the harmonic oscillator. 

IV. CONCLUSIONS 

The results of the previous sections confirm the findings of references [§,0, and gener- 
alize them to arbitrary D. The semiclassical approach finds the minima of the Euclidean 
action and expands around them. As a result, it generates a series whose terms correspond 
to resummations of infinite numbers of perturbative graphs plus additional ones. Our cal- 
culations show that even the lowest order semiclassical estimates improve on perturbation 
theory at low temperatures and, in contrast with it, correctly describe the high temperature 
regime. 

The comparison with WKB estimates, done for the one- dimensional case, is particularly 
interesting. Such estimates approximate the values of the energy levels of the single-well 
anharmonic oscillator to a high precision ii g = h\/m?uj^ is small, even if we restrict ourselves 
to the lowest order WKB quantization condition, given by the Bohr-Sommerfeld formula 
They were then used to compute the partition function by actually performing the 



sums over eigenstates numerically. Thus, the WKB results can be considered "quasi-exact" . 
By contrast, the semiclassical approach directly approximates the whole sum. Its lowest 
order agrees well with the quasi-exact WKB result at both high and low temperatures. At 
high T, this agreement just reflects the convergence of both results to the classical limit, 
something which is completely missed by perturbation theory. Only in the intermediate 
region does our result differ from WKB, although we expect this to be modified with the 
inclusion of next-to-leading orders. It is less accurate, as it approximates the whole sum, 
whereas WKB approximates each term in the sum; but it does incorporate and improve 
upon the virtues of perturbation theory at low temperatures, and of the classical limit at 
high ones. Results for D = 2 and D = 3 do follow the same pattern, although we have not 
compared them to WKB estimates. 

The advantage of this method is that it reduces the whole quantum problem to the 
computation of (few) classical paths. From then on, a systematic procedure takes care of 
generating each term in the series. Paradoxically, this may also be its weakness: there are 
systems for which the action does not have a global minimum, but which are perfectly well- 
defined quantum-mechanically. The Coulomb potential is a good example; there, depending 

9 



on the values of {3 and tq, the number of classical paths may be two, one or zero. Besides, 
only in the two-solution regime do we have minima. Even then, they are local, not global 
ones. Therefore, our starting point seems ill-defined. This should not come as a surprise, 
however, since here the classical limit itself is ill-defined, as the potential is unbounded 
below. As a matter of fact, even the usual time-slicing prescription to calculate the path 
integral must be modified in the case of the Coulomb potential [§]. Cases like this will require 
special consideration, although there exist suggestions in the literature as to how to treat 



similar situations of absence of classical paths in Quantum Mechanics |T9| . Nevertheless, we 
expect the techniques presented here to be useful in any problem which can be reduced to 
the calculation of partition or correlation functions in equilibrium Statistical Mechanics, as 
long as it allows for a simple analysis of the minima of the Euclidean action. 

Our next step is to investigate how the semiclassical treatment affects field-theoretic 
problems at finite temperature, where standard methods of computation of effective poten- 
tials rely on expansions around constant backgrounds. At finite temperature, these are not 
in general minima of the Euclidean action. This might lead to problems with the expansions 
around such backgrounds at high temperatures, of the same nature of those encountered by 
perturbation theory in Quantum Statistical Mechanics. Even if we neglect any coordinate 
dependence of the fields, their dependence on Euclidean time is essential to satisfy equations 
of motion and boundary conditions that caracterize classical paths. We expect this to have 
an effect on a variety of calculations. 

Another problem of interest is to generalize our results to field theories with spherically 
symmetric classical solutions. An extension of the approach presented in this work to treat 
models containing non-trivial backgrounds (like instantons, monopoles, vortices, etc.) as 
classical solutions might lead to some new insights. Unfortunately, the extension of our 
results to field theories is not a straightforward process. In fact, we do not know how to 
construct a semiclassical propagator in general. The success of our program will depend on 
how well can we circumvent this difficulty. 
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APPENDIX A: 

Let J{t, t') be the solution of the homogeneous differential equation 






satisfying the initial conditions 

10 



J,-,(r,r') = (Al) 



^(r',r') = 0, l-J(r = T',T' 



1 
— 1 
m 



(A2) 



This function is known as the Jacobi commutator ||TT],|T2|. It can be exphcitly constructed 
as follows. Let x(r; a, b) be the solution of the equation of motion (0) satisfying the initial 
conditions x(0) = a, x(0) = b. Let A and B be the D x D matrices defined as 

d 

Ajkir) = -—Xj{t; a = Xo,h = Vo), (A3) 



dttk 



d 



^jk{r) = — - Xj{t; a = Xo, b = vq) 



(A4) 



where vg = Xc(0). By differentiating Eq. (^ with respect to a^ and bj. (and taking a = Xq, 
b = Vq), one can show that they are solutions of ([Al|). They are also invertible for r 
small enough |^ (but not zero). Indeed, x(r) = a + br + O^r"^) when r — » 0, hence 
A{t) = 11 + O(t^) and B{t) = r U + O(r^). Therefore, the expression 



J{t, r') = -— \A{t) A'\t') - B{t) B-\t')] \A{r') A'^r') - B{t') B-\t')' ' 

ffl I J L 



(A5) 

makes sense, and one can easily verify that it satisfies ( [Al| ) and ([A2|) . 

The Green's function G{t, t') can be written in terms of the Jacobi commutator as 

G{t, t) = J{t, 0) M(0, (3n) J{(3h, t) 6{t' - t) - J{t, (3h) M{(3h, 0) J(0, r ) 0(r - r') , (A6) 

where M{t,t') = —J{t',t)^^ and 6{t) is the Heaviside step function. To prove (|A6|) we 
need the following identities: 



J{t, 0) M(0, (5h) J{f3h, t') + J{t, ph) M{f3h, 0) J(0, r') = -J{t, t') 

1 



drJir, 0) M(0, (3h) J{(3h, r) + d^Jir, (3h) M{(3h, 0) J(0, r) 



1 . 



TTl 



(A7) 



(A8) 



The first identity follows from the fact that both functions are solutions of the same second 
order differential equation [Eq. ([AID] and are equal at r = and r = j3h; the second follows 
from dAD and (^). 



Now, the proof of (|A6|) : (i) it is a solution of (^ when r < r' or r > r'; (ii) it satisfies 
the boundary conditions ([10|); (iii) it is continuous at r = r', 

^(r' + 0, r') = G(r' - 0, r') (A9) 

[use (|A7|) with t = t'], and (iv) its derivative with respect to r has the discontinuity implied 
by(D, 



I- G(r = r' + 0, r') - A G(r = r' - 0, r') = --i 

OT OT m 

[use (^]. 

Finally, the determinant A of the fluctuation operator JF is given by 

A = (27r;i)^ det [- J(/?^, 0)] . 

This formula can be proven along the lines of Appendix 1 of Ref. [pT . 



(AlO) 



(All) 
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FIGURES 

FIG. 1. Specific heat (in units of ks) vs. temperature (T = 1/ (3Tiuj) for the one- (diamonds), 
two- (circles) , and three-dimensional (crosses) single- well quartic anharmonic oscillator in the semi- 
classical approximation, g = h\/m?uj^ = 0.5. 

FIG. 2. Specific heat (in units of ks) vs. temperature (T = l/f3hcv) for the one-dimensional har- 
monic oscillator (long-dashed line) and for the single-well quartic anharmonic oscillator: classical 
result (short-dashed line), semiclassical approximation (circles), and WKB approximation (solid 
line), g = ?i\/m?uj'^ = 0.2. 
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